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We combine recent measurements of Cosmic Microwave Background Anisotropics, Supernovae 
luminosity distances and Baryonic Acoustic Oscillations to derive constraints on the dark energy 
equation of state w in the redshift range < z < 2, using a principal components approach. We find 
no significant deviations from the expectations of a cosmological constant. However, combining the 
datasets we find slight indication for w < — 1 at low redshift, thus highlighting how these datasets 
prefer a non-constant w. Nevertheless the cosmological constant is still in agreement with these 
observations, while we find that some classes of alternative models may be in tension with the 
inferred w(z) behaviour. 

PACS numbers: 98.80.Es, 98.80.Jk, 95.30.Sf 



I. INTRODUCTION 

One of the main goal of modern cosmology is to de- 
termine the nature of the dark energy component that is 
sourcing the late time accelerated expansion of the uni- 
verse. 

Since its discovery from measurements of luminosity 
distances of type la Supernovae in 1998 [l|, H|, cosmic 
acceleration has been confirmed by several independent 
cosmological data. In particular measurements of the 
Cosmic Microwave Background and galaxy distribution, 
have not only confirmed its presence but also helped in 
clarifying its nature. 

Constraints on the dark energy equation of state w, 
that is the ratio between dark energy pressure and den- 
sity, assumed as constant with redshift, have slightly im- 
proved over the last decade. For example, from the early 
analysis of [|| where the constraint of w < —0.6 at 95% 
c.l. was reported, the very recent analysis of Q gives 
w = —1.01 ± 0.085, at 68% c.l., i.e. an improvement 
of a factor 2.35 in the constraint in about a decade of 
observations. These recent observations have provided 
decisive evidence against dark energy models with values 
of w in the range —0.8 < w < —0.66 as domain walls (see 
ej;. Q) or models based on extra dimensions (see e.g. 
Q and references therein). Other models as "geometric" 
dark energy and "thawing" dark energy show, albeit not 
at a statistically significant level, a better agreement with 
observations (see Q). 

Up to now, data are in good agreement with the 
standard cosmological model, where the acceleration 
is sourced by the cosmological constant A Q, and 
therefore the equation of state is w = —1. However, 
recent experiments have reached a sensitivity that allows 
to test other characteristics of the dark energy equation 
of state and not only its value at very recent times. It is 
indeed possible to reconstruct the redshift dependence of 
w, a key feature to distinguish between the cosmological 
constant, which predicts a constant w, and alternative 
models which generally give an effecitve equation of state 



parameter that varies with redshift (see e.g. [3 El). 

In this paper we investigate the power of currently 
available cosmological data to constrain w(z) at different 
redshifts, adopting the principal component analysis 
already applied to the equation of state of dark energy 
in previous studies (see e.g. (SrtHj]). 

The paper is organized as follows. In Section [TTJ we 
list the datasets used in our analysis, while in Section IILTI 
we describe the methodology. The results of our Markov 
Chain - Monte Carlo analysis are presented in Section ITVl 
while we draw our conclusions in Section [V] 



II. DATA 

We analyze a large sample of cosmological datasets. 
For the supernovae SN-Ia luminosity distance we con- 
sider the Unionl compilation [15], Union2 [l6j and 
SNLS [13, Ell ■ We then combine them separately with 
the CMB observations coming from seven years of obser- 
vations from the WMAP satellite [|]. 

Together with the WMAP and SNLS we consider dif- 
ferent datasets for Baryonic Acoustic Oscillations (BAO) 
surveys combined in the following way: 

• runl 

SDSS-dr7 at z=0.20, 0.35 [11 in form of d s (z), Wig- 
gleZ at z=0.44, 0.60, 0.73 [2Cj in form of A s (z) ; 

• run2 

6dFGRS at z=0.1 [Uj], WiggleZ at z=0.44, 0.60, 
0.73 [H all in form d s (z); 

• run3 

WiggleZ at z=0.44, 0.60, 0.73 all in form d s {z); 

• run4 

6dFGRS at z=Q.l, SDSS-dr7 at z=0.20, 0.35, Wig- 
gleZ at z=0.44, 0.60, 0.73 all in form d s (z). 
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In runl we use as covariance matrix the one indicated 
in for SDSS-dr7 points and the one shown in [2(| for 
WiggleZ points. In the other runs we use always the co- 
variance matrix of 

.19] 

for SDSS-dr7 where used, while for 
6dFGRS and WiggleZ we use their parts of the covariance 
matrix from p3| . Finally we add the most recent mea- 
surements for Hq from the Hubble Space Telescope 24 1 
and the H(z) dataset from Moresco et al. [25j . 



III. DATA ANALYSIS METHOD 

The analysis method we adopt is based on the publicly 
available Monte Carlo Markov Chain package cosmomc 
|26j with a convergence diagnostic done through the Gel- 
man and Rubin statistic. 

We sample the following six-dimensional standard set 
of cosmological parameters, adopting flat priors on them: 
the baryon and cold dark matter densities fib and f2 c , the 
Hubble constant Hq, the reionization optical depth r, the 
scalar spectral index n<j, and the overall normalization 
of the spectrum As at k = 0.002 Mpc~ . We consider 
purely adiabatic initial conditions and we impose spatial 
flatness. 

The dark energy equation of state, as discussed above, 
is sampled in six redshift bins, Wi(z) (i — 1,2,. .6), at 
six redshifts, z, e [0.0,0.25,0.50,0.85,1.25,2.0] equally 
spaced in ln(a). Including more than six bins does not 
significantly improve the constraints. 

In order to have w(z) as a smooth and continuous 
function, we interpolate between these values with a 
hyperbolic tangent function defined as: 

w(z) — Wi Z = Zi 

w(z) = Wi + S w + S w tanh(^f^) z e [z t , z i+1 ] 

w(z) = — 1 Z > Z$ 

where Ai is an amplitude factor, S z is the half value of 
(zi + i — Zi), S w is the half value of (w;+i — Wi) and s is a 
smoothing parameter. The presence of dark energy per- 
turbations has to be accounted for when the dark energy 
equation of state is constrained using cosmological probes 
that are sensitive to density perturbations. This is done 
by using a modified version of the publicly available code 
CAMB [29|, that evolves the dark energy perturbations 
also for values of w < —1, avoiding singularities by using 
the Parametrized Post-Friedmann prescription for dark 
energy suggested by Hu [3(il - [32l |. 

Once the values of the parameters w = Wi are deter- 
mined, we must deal with the fact that they are cor- 
related, which means that their covariance matrix is 
not diagonal. To obtain an uncorrelated estimation for 
their value we can make a rotation in the parameters 
space so to chose a different basis for the w in which 
their covariance matrix is diagonal. This can be done 
following the procedure suggested in [{J [Io|, l27j . Us- 
ing CosmoMC we derive the covariance matrix, C = 
(Wi — {wi))(wj — {wj)) T = (ww T ) — (w)(w T ), and then 
we invert it to obtain the Fisher matrix, which can be 



rewritten also as C _1 = F = O t AO. Here A is the 
diagonalized inverse covariance for the transformed bins. 
The vector of the uncorrelated parameters q is obtained 
using the rows of the tranformation matrix as weights 
q = Ow. We can now define W so that W T W = F, 
and, as pointed out by Hamilton [28] , there are infinitely 
many choices for W. Following Q we choose as weight 
matrix W = T A2 0, where the rows are normalized to 
unity, and we apply it to obtain the uncorrelated param- 
eters. 

In other words we run CosmoMC with the w vector 
to obtain its covariance matrix. After that we derive 
q, the principal components vector, by using the proce- 
dure explained before, and then run again the statistics 
to evaluate its marginalized values and errors. In Figure 
[I] we show typical weights to obtain the q from a linear 
combination of w while in Figure [5] we show the like- 
lihoods of the six principal components qi for a typical 
run. It can be noticed how, as expected, the lower red- 
shift components are better constrained than the higher 
redshift ones. 
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FIG. 1. Typical weights used to obtain the principal com- 
ponents qi starting from the values of the w%. The weights 
plotted are those obtained by PCA decorrelation for run2 of 
Table HV1 



IV. RESULTS 
A. WMAP+SNIa 

We first investigate the constraints on the dark energy 
equation of state combining WMAP with different super- 
novae datasets and we report the results of this analysis 
in Table |H] and Figure U 

As we can see, all the values for the qi are consistent 
with the cosmological constant case with qi = — 1 in be- 
tween two standard deviations. However, if we look at 
the value of 53, around z ~ 0.5, the SNLS survey ap- 
pears as more compatible with w = — 1 than the Unionl 
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FIG. 2. Typical likelihoods for the principal components qi. 
This results are obtained with PCA decorrelation for run2 of 
Table IW1 
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FIG. 3. Constraints at 95% c.l. on q t from WMAP CMB data 
combined with different supernovae catalogs: Unionl, Union2 
and SNLS. As we can see all the datasets provide constraints 
on qi that are consistent with the predictions of a cosmological 
constant in between two standard deviations. 



and Union2 catalogs that are on the contrary provid- 
ing slightly larger values of q%. Moreover the SNLS 
gives stronger constraints at lower redshifts respect to 
Unionl and Union2 while appears preferring a value of 
q4 < — 1 at z ~ 0.85 whereas Unionl and Union2 are 
more compatible with a cosmological constant in this red- 
shift range. 

If we look at the constraints on the other parame- 
ters in Table [TT] we also see that the three datasets pro- 
vide slightly different values for the Hubble constant Hq: 
the SNLS survey indicates larger values, with Hq ~ 
72.5km/s/Mpc, while the Union2 and Unionl catalogs 
prefers smaller values with Hq ~ 70Akm/s/Mpc and 
Ho ~ 68.2km/s/Mpc respectively. All these values of 
Hq are however well consistent in between two standard 
deviations. The value of the dark energy density appears 
also as larger in the SNLS survey respect to the value ob- 
tained from the Union2 and Unionl datasets, the latter 
providing the smallest value. Since we are considering a 
flat universe this means that SNLS data is preferring a 
lower matter density respect to Union2 and Unionl. 



B. WMAP+SNLS+BAO 
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FIG. 4. Constraints on the dark energy equation of state from 
a combined analysis of the WMAP+SNLS dataset with differ- 
ent BAO datasets as described in the text. The run2 provides 
the lowest chi-square value. There is a broad agreement with 
a cosmological constant, however values at redshift z ~ 0.25 
are in better agreement with w < — 1 while values at z ~ 0.5 
prefer w > — 1. 



Since the SNLS dataset appears as the most consistent 
with a cosmological constant we perform several analysis 
combining the WMAP+SNLS data with four different 
choices for the BAO datasets as specified in the previ- 
ous Section. The results on the parameters are reported 
in Table IIVI while the constraints on qi are reported in 
Figure |U 

We can derive the following conclusions: 

• There is a slight indication for w < — 1 around 
the second bin q 2 at z ~ 0.25. This deviation, 
already noticed in [l2j], is at about two standard 



deviations and appears in all the BAO data combi- 
nations considered. The values of qs in the adjacent 
bin at z ~ 0.50 are on the contrary larger than — 1 
at about one standard deviation. The agreement 
with a cosmological constant is worse respect to 
the WMAP+SNLS case. 

• When the BAO data is included, the preferred 
value for the Hubble constant is significantly 
smaller respect to the WMAP7+SNLS case. The 
values for Hq are in the range 64 — 65km/s/Alpc. 
These values are in strong tension with the HST 
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TABLE I. WMAP+SNIa 



T* E» r E» TY1 or 




TTnirm9 






0.0222 ±0.0011 


0.0222 ±0.0011 


0.0222 ± 0.0011 


Q c h 2 


0.113 ±0.011 


0.113 ± 0.011 


0.113 ±0.010 


Ho 


68.2 ±6.1 


70.4 ±6.1 


72.6 ± 7.0 


n s 


0.963 ±0.028 


0.963 ± 0.026 


0.962 ± 0.027 


log[W 10 A a ] 


3.209 ±0.091 


3.212 ± 0.085 


3.213 ± 0.084 


Asz 


1.0 ± 1.2 


1.0 ± 1.2 


0.9 ± 1.1 


a,-(z — 2 OO 1 ) 


> —3.0 


-1 6+H 


— 1.5 ± 1.3 


nz.(r — 1 9^ 
f/5 1,2 — 

qi{z = 0.85) 
g 3 (z = 0.50) 
92(2 = 0.25) 
qi{z = 0.00) 


1 - 6 -1.2 

-0.92+^2 
-0.84+^° 
-1 02+ 022 

n OO+0.39 


-1.11+"" 

— U. 71 

-0.89±° -Ji 
-1 02+ 029 

1-U^_ 28 


1 ai+0.80 
— 0.87 

-1.06 + °- 3 37 

-1.03 ±0.21 


n A 


0.708 ±0.061 


0.725 ± 0.059 


0.742 ± 0.057 




13.84 ±0.32 


13.79 ±0.31 


13.72 ±0.32 




0.292 ±0.061 


0.275 ± 0.059 


0.258 ± 0.057 




1.0380 ± 0.0054 


1.0380 ± 0.0052 


1.0382 ± 0.0052 




7776.58 


7999.48 


7586.84 



TABLE II. Constraints at 95% confidence level for a WMAP analysis considering different supernovae datasets. The SNLS 
survey provides constraints that are more consistent with the A-CDM case. 



TABLE III. WMAP+SNLS+BAO 



Parameter 


runl 


run2 


run3 


run4 


Q c h 2 
Ho 

n s 

log[10 10 A s ] 

A S z 
q 6 (z = 2.00) 
qs(z = 1.25) 
q4z = 0.85) 
93(2 = 0.50) 
q 2 (z = 0.25) 
qi(z = 0.00) 


0.0220 ±0.0011 
0.119 ± 0.012 

64.3 ±3.5 
0.957 ± 0.029 
3.245 ± 0.098 
0.9 ±1.1 
> -2.7 

-0-841 

i O1 +0.l8 
L - zl -0.19 

-1 06 +0 - 27 

i - uu -0.30 


0.0221 ± 0.0011 
0.116 ±0.011 

65.5 ±3.8 
0.960 ± 0.028 
3.229 ± 0.093 
1.0 ± 1.1 
-1 10+ 098 

-0.86+HI 
-0 fi4+° :4 ° 

u - D *-0.68 

-0.80lg;f 9 

i 1O+0.18 
^^^-O.M 
-, n o+0.27 


0.0221 ±0.0011 
0.117 ±0.012 

64.4 ±4.2 
0.960 ± 0.029 
3.229 ±0.095 
1.0 ± 1.1 

-0.67±°;1| 

n 47 +0.34 
u - 4 '-0.62 

-0 «^+°- 28 

_1 14+018 
-1.08+0 29 


0.0222 ± 0.0012 
0.117 ±0.012 

63.6 ±3.0 
0.962 ± 0.031 
3.22 ±0.10 
1.0 ± 1.1 

-0.79+?^ 
-0.59+^ 
-0-45±C;^ 

n 7 o+0.26 
u - lo -0.37 
i i 7 +0.18 
L - Ll -0.20 

-1 06+ ' 27 

1 - UO -0.32 


n A 

to/Gyr 

e 


0.659 ± 0.035 

14.1 ±0.28 
0.341 ± 0.035 
1.0370 ± 0.0052 


0.676 ± 0.040 
14.01 ±0.28 
0.324 ± 0.040 
1.0375 ± 0.0054 


0.665 ± 0.048 
14.08 ± 0.30 
0.335 ± 0.048 
1.0373 ± 0.0054 


0.656 ± 0.035 
14.11 ±0.28 
0.344 ± 0.035 
1.0374 ± 0.0055 


x 2 


7595.08 


7591.30 


7593.98 


7597.02 



TABLE IV. Constraints at 95% confidence level for a WMAP±SNLS analysis considering different combinations of BAO 
datasets (see text). BAO data clearly prefer a lower Hubble constant around Ho ~ 65km/s/Mpc. 



result of Hq = 73.8 ± 2Akm/s/Mpc of [24j and • The BAO datasets are all giving consistent results, 

may indicate the need for new physics or, more However the "run2" case provides the lowest best 

simply, the existence of systematics. fit chi-square. 



• Similarly, the value of the matter density is 
larger when the BAO dataset is included in the 
WMAP±SNLS dataset. 
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C. WMAP+SNLS+BAO+HST and +H(z) 

Since the run2 BAO dataset provides the low- 
est chi-square values, we assume it, together with 
WMAP+SNLS, as our basic dataset and we now add 
the HST prior on the Hubble parameter and the H(z) 
dataset. The results are reported in Table I VII and in 
Figure [5] 
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FIG. 5. Constraints on the dark energy equation of state from 
a combined analysis of the WMAP+SNLS+BAO dataset with 
the HST prior on the Hubble constant of [24| and the determi- 
nation of H(z) in [25j]. There is an improved agreement with 
a cosmological constant respect to the WMAP+SNLS+BAO 
case for 53 at z ~ 0.50. 

As we can see the effect of adding the HST prior or 
the H(z) dataset is to increase the value of the Hubble 
parameter. However the HST prior is clearly in tension 
with the WMAP+SNLS+BAO dataset as showed by the 
increased value of the chi-square. 

Both the HST prior and the H(z) dataset render the 
value of q% more compatible with predictions of a cosmo- 
logical constant (w = —1). In general, with the exception 
of the value of q<± at z ~ 0.25 that prefers values such that 
W < — 1, there is a general agreement with a cosmological 
constant. 



D. Comparison with alternatives to ACDM 



our results which slightly prefer w < — 1; the covariant 
galileon model instead does give w < — 1 however, while 
it can fit our first two redshift bins, it has tension with 
values in the higher redshift bins. 

We want to stress that in this paper we did not obtain 
constraints on these theories; our purpose is to show how 
a detailed analysis of currently available datasets can in 
principle rule out models that are still considered viable. 
In particular, while we know that the covariant galileon 
model is already disfavoured by other observations (36[, 
tracking models and HS modified gravity are still viable 
for some interval of the parameter space [13, HU . How- 
ever, they appear to be in tension with the combination 
of datasets considered in our analysis, motivating a de- 
tailed full analysis that will be the matter of an upcoming 
paper. 
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FIG. 6. Comparison of the constraints on the dark energy 
equation of state derived from a WMAP+SNLS+BAO+H(z) 
dataset with several dark energy models as covariant galileon, 
Hu-Sawicki, tracking SUGRA and tracking power-law. Most 
of these models are in tension with the data mostly because 
of the low redshift values of q\ and q2 that point to a cos- 
mological constant or to w(z ~ 0.25) < —1. The tracking 
power law model reproduces the low redshift behaviour but 
is in disagreement with the constraints at higher redshifts. 



In order to better picture the impact of our analysis we 
compare the constraints on the equation of state obtained 
using WMAP+SNLS+BAO+iJ(z) reported in the pre- 
vious Section, with some representative dark energy and 
modified gravity models alternative to ACDM. 
In FigJH] we show how our results on qi are in tension 
with the typical w{z) behaviour of models such as the 
Hu-Sawicki f(R) model (HS) [30|, the covariant galileon 
model [H[ and tracking models|34l |35| . It is possible to 
notice how the tension is related to the low redshift bins; 
in particular, we can see how both tracking and HS mod- 
els predict a w(z) > — 1 at low redshifts colliding with 



V. CONCLUSIONS 

In this paper we have presented updated constraints on 
the dark energy equation of state from a large collection 
of cosmological datasets. We have found that while for 
a WMAP+SNLS case a cosmological constant is in very 
good agreement with the data, when also BAO data are 
considered some indications are present for a time evolu- 
tion of the dark energy equation of state. In particular 
the values around z ~ 0.25 prefer w < — 1 above one 
standard deviation. 
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TABLE V. WMAP+SNLS+BA0(run2) 



Parameter 


+HST 


+H(z) 


±H(z)±HST 




0.0222 ±0.0011 


0.0219 ±0.0011 


0.0221 ± 0.0011 


Q c h 2 


0.121 ±0.012 


0.1248 ±0.0096 


0.1250 ±0.0093 


Ho 


67.2 ±3.4 


66.4 ±3.0 


68.1 ± 2.7 


n s 


0.959 ±0.029 


0.950 ± 0.027 


0.953 ± 0.026 




9 O/IQ _(_ n riQ7 
O.Z^io it u.uy / 


^ 98/1 -4- n 081 


9 977 -1- fl 077 


Asz 


i n i 1 9 
l.U rt X.z 


n n 1 i i 


n q 4- i i 
U.y ± l.l 


qe(z = 2.00) 
q 5 {z = 1.25) 

n ,( — n 

g 3 (^ = 0.50) 
q 2 (z = 0.25) 
gi(z = 0.00) 


> —2.9 

n 7n+ - 56 
-1.04+21 

— (J.4» 

_] iq+0-19 

_1 1 4+0-29 


1 9 + 1-1 
1 O1+0-87 

-l- 21 -i.ia 

i m+°- 53 
-i.oi+S-22 

— U.ob 

1-1O-0 .18 
-i -19+0.23 


1 9 + 1-1 

-J--2-1.5 
1 Qfi+0-95 
l.oD_ 9g 

1 is+o- 57 

1-18_0.69 

-1.05+S-S 

— (J.o4 

-1.18 ±0.17 
-1 16+ 025 


n A 


0.683 ± 0.036 


0.667 ± 0.037 


0.682 ± 0.033 


ta/Gyr 


13.92 ±0.24 


13.96 ± 0.24 


13.86 ±0.22 




0.317 ±0.036 


0.333 ± 0.037 


0.318 ±0.033 




1.0380 ± 0.0052 


1.0368 ± 0.0054 


1.0376 ± 0.0052 


x 2 


7608.70 


7611.00 


7619.64 



TABLE VI. Constraints at 95% confidence level for a WMAP+SNLS+BAO analysis considering a HST prior on the Hubble 
constant of [24| ] and the determination of H(z) in 25] 



Moreover, BAO data, when analyzed in our extended 
parametrization of the dark energy component, pre- 
fer a low value for the Hubble constant with Hq ~ 
65km/s/Mpc in clear tension with the measured HST 
value of 

Finally, we have compared our results with the theo- 
retical predictions of dark energy and modified gravity 
models that predict a time varying equation of state, 
finding that several models do not reproduce the w(z) 
trend obtained with our analysis. This comparison sug- 
gests that the combination of datasets considered in our 
analysis can offer powerful constraints on alternatives to 
the cosmological standard model. We intend to perform 



a full analysis of the reduction of the allowed parameters 
space for the latter, when they are constrained directly, 
in an upcoming paper. 
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